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1. Introduction 

In this contribution we present an analysis of the static potential on lattices generated with 
Wilson plaquette gauge action and Nf = 2 flavors of non-perturbatively 6 {a) improved Wilson 
quarks. Periodic boundary conditions are used for all fields apart from anti-periodic boundary 
conditions for the fermions in time. The static potential is extracted from measurements of the 
on-axis Wilson loops and can therefore only be determined for distances smaller than the string 
breaking distance. 

2. Method 

Our method to extract the static potential from HYP smeared Wilson loops has been presented 
in §\§ and is based on two steps. In the first step all the gauge links are HYP smeared [@] using the 
HYP2 parameters 

ai = 1.0, a 2 = l-0, a 3 =0.5. (2.1) 

As the time-like links are concerned, this is equivalent to the choice of a static quark action. The 
binding energy of a static-light meson is Z? stat ~ + • • • an d the HYP2 choice of smearing 

parameters minimizes e^ 1 ' [Q]. The static energies V n (r) (Vb = V is the static potential) depends on 
the static quark action. Their contribution to the expectation values of Wilson loops W(r, T) of size 
r in one of the spatial direction and size T in the time-direction is given by 

(W(r, T)) ~ £c„c*e- y " W(r - 2fl) (with N t -> °o time-slices) . (2.2) 

n 

The coefficients c„ depend on the choice of the space-like links. In the second step of our procedure 
we construct a variational basis for them. The space-like links are smeared using ni iterations of 
spatial HYP smearing 

(W(r,t)) — > C lm (r,T) (2.3) 

We use a basis of I = 1,2,3 operators and the generalized eigenvalue method to extract V n [Q, [5]]. 
The error analysis is based on the method of [||] and we add a tail to the autocorrelation function to 
account for the slow modes 

In this contribution we mainly present results for two lattice sets at approximately constant 
pseudo-scalar mass ryMps pa 0.62 . . . 0.64: 

F7a: 701 configurations , 96 x 48 3 , r /a(F7a) = 7.079(63) , (2.4) 
07: 740 configurations , 128 x 64 3 , r /a(O7) = 9.63(12) . (2.5) 

In the left plot of Fig. [I] we show the two potentials renormalized by subtracting twice the energy 
of a static-light meson. The right plot of Fig. [j] compares the potentials on the A5a, E5g and N5 
lattices [JsJ] computed with two different sets of HYP parameters for the static action: the HYP2 set 
in Eq. ( |Q| ) and the set called HYP [g]. The difference of the potentials is 

ro[V HYP2 (r)-y HYP (r)] = 2r (*™ -E™) + ^G(Ar,^r) . (2.6) 
The dotted lines and gray bands represent the values with errors of 2 ro {E^^ 2 — E^"). We observe 



small cut-off effects encoded in the function G in Eq. ( |2.6| ) 
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Figure 1: Left plot: the static potentials on the F7a and 07 lattices. Right plot: a comparison of the potentials 
using two static actions on the A5a, E5g and N5 lattices. 

3. Scale ro/a 

Here we present a preliminary update on the determination of the scale r [§] at the three 
available j8 values and based on the currently available statistics. Motivated by the right plot of fig. 
H] we perform a combined linear extrapolation in the square of the pseudo-scalar mass to the chiral 
limit and obtain the following results for ro/a(j8): 

r /a(5.2) =6.159(78), r /a(5.3) = 7.242(70) , r /a(5.5) = 9.84(13) . (3.1) 

With these values of ro/a the preliminary update of the ALPHA result for the A parameter extracted 
from the running of the Schrodinger Functional coupling jlQ| ] is 



0.78(3)(5). (3.2) 

The final numbers will be given in [Hi. Alternative methods for determining the lattice spacing in 
the CLS consortium are presented at this conference, through the kaon decay constant fk [11] or 
the Q. mass [|2|]. 

4. Mass dependence and cutoff effects in ro 

In order to discuss the mass dependence and cut-off effects in ro, we introduce the variable 
x = (roMps) 2 , where Mps is the pseudo-scalar mass, and define a reference value ro| re f, which 
corresponds to the value of ro at the pseudo-scalar mass 

* re f = (roMps) 2 | re f = 0.75. (4.1) 
In the left plot of Fig. ^| we plot ro / ro | re f versus x and compare our data with the twisted mass data 



of [|L3j. The value of ro| re f is obtained by linear interpolation and the error analysis of ro/ro| re f 
takes into account the correlations between the data. We only consider data with x < 1.2. 

Using the data plotted in the left plot of Fig. we determine for each j3 value the slope 
s(a/ro| re f) in the Taylor expansion 



ro ref 



-(x) = l+j(a/ro|ref)-(x-JCref). 



(4.2) 
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Left plot: mass and cut-off dependence of ro. Right plot: continuum extrapolation of the slope 
Comparison with twisted mass data from [JljJ]. 



The results are shown in the right plot of Fig. ||| In principle they depend on the chosen value of 
*ref = 0.75. We repeated the analysis using x ve f = 0.6 but the slopes do not change significantly. 
The twisted mass results for s(a/ro| re f) increase in magnitude as a — > and they line up with our 
results at the smallest a. Our data do not show significant cut-off effects and we emphasize that 
the errors are larger due to accounting for the slow modes in the autocorrelation function. In the 
continuum limit we obtain the value s(0) = —0.0658(95) by fitting to a constant. 

5. Running couplings 

The static force F(r) = V'(r) defines the running coupling in the qq-scheme g qq through 



2 ( At ) = ^r 2 F(r), M 



1/r. 



(5.1) 



We use an improved lattice definition of F which removes the cut-off effect effects at tree level 
[jl|]. The results for a qq (l/r) = g^ q / (An) on the F7a and 07 lattices are plotted in the left plot of 
Fig. ||. The magenta dashed and continued lines are the 3- and 4-loop running respectively. They 
are obtained by solving the renormalization group (RG) equation, which for a given scheme S reads 
(here S=qq) 



As 



exp 



1 



+ 



1 



j3 s (x) box 



(5.2) 



For the beta function /3s we use the 3- and 4-loop approximations jl8| , [P^ , [20| , and roAs is 
calculated from the value in Eq. ( p\2[ ) [p. One sees that our smallest lattice spacing is not small 
enough to make contact with perturbation theory. This becomes more evident by the attempt to 
determine the A parameter through Eq. (5.2) shown in Fig. ^ We use the a qq (l/r) values measured 
on the 07 lattices as input for the right hand side of Eq. ( |5.2| ) together with the perturbative beta 
function (2-loop, filled circles; 3-loop filled squares; 4-loop filled diamonds) and calculate foA^. 
For a qq (1/r) — y the values for the different loop approximations should converge to one value, but 
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Figure 3: Left plot: the coupling a qq (l/r) compared to perturbation theory and Nf = results from QljJ], 
Right plot: the slope c(r) compared to perturbation theory, potential models and iVf = results from ||15||. 



this is not the case for our data in Fig. |] (note that there is no sign of the 2-loop data bending toward 
the 3- and 4-loop ones). For comparison we plot the Nf = data from jl4| ] which demonstrate that 
couplings down to at least a qq (l/r) pa 0.2 (half the value of our smallest coupling on the 07 lattice) 
would be needed. 
From the slope 



c( r ) = -r 3 F\r) 



the running coupling g c is obtained 



'(AO 



4n 



c(r), ^ = 1/r, 



(5.3) 



(5.4) 



which defines the c-scheme. The results on the F7a and 07 lattices are plotted in the right plot 
of Fig. [| We also plot the Nf = data from [15] showing that there is a large effect originating 
from dynamical fermions (at distances much below the string breaking distance), as already indi- 
cated by the perturbative curves. For comparison we plot the values of c{r) calculated using two 
phenomenological potential models: the Cornell potential [ |l6| ] 



Vci 

and the Richardson potential [ |l7j ] 

Vk(r) 



Var, fc = 0.52 

r 



1 -A ( Ar f(Ar) 



6nbo 



Ar 



(5.5) 



(5.6) 



If rA 1, Vfi(r) ~ —l/[6nbor\n(l/Ar)], where bo is the 1-loop coefficient of the beta func- 
tion. If rA ^> 1, Vr(t) ~ const x r. Our data seem to follow the Richardson curve (dotted line, 
Nf = 2, roA = 0.73) and approach the Cornell value at distances around ro, where unfortunately the 
statistical errors (enhanced by the factor r 3 ) become too large. 
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Figure 4: Attempt to extract the A parameter solving the RG equation for a qq (l/r) with the 2-, 3- and 
4-loop /3 function. Nf — data are from | 



6. Conclusions 

We determine the static potential from Wilson loops using the HYP2 static action. Cut-off 
effects in the potential appear to be small and ro/a can be accurately determined. There could be 



non-negligible cut-off effects in the mass dependence of ro [22], however they seem to be small 
with improved Wilson fermions. 

The running coupling from the static force is compared to perturbation theory: at our lattice 
spacing we cannot reach small enough distances in order to extract the A parameter. This is in 



contrast to what is concluded by a recent work with the same lattice spacing [J23[] . We observe large 
effects from dynamical fermions in the slope c{r). The shape of the static potential is relevant for 



holographic QCD models [24, 25]. 
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